Signatures of Hong-Ou-Mandel Interference at Microwave Frequencies 
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Two-photon quantum interference at a beam splitter, commonly known as Hong-Ou-Mandel in- 
terference, was recently demonstrated with microwave-frequency photons by Lang et al. [T . This 
experiment employed circuit QED systems as sources of microwave photons, and was based on 
the measurement of second-order cross-correlation and auto-correlation functions of the microwave 
fields at the outputs of the beam splitter. Here we present the calculation of these correlation func- 
tions for the cases of inputs corresponding to: (i) trains of pulsed Gaussian or Lorentzian single 
microwave photons, and (ii) resonant fluorescent microwave fields from continuously-driven circuit 
QED systems. The calculations include the effects of the finite bandwidth of the detection scheme. 
In both cases, the signature of two-photon quantum interference is a suppression of the second-order 
cross-correlation function for small delays. The experiment described in Ref. [T|was performed with 
trains of Lorentzian single photons, and very good agreement between the calculations and the 
experimental data was obtained. 



I. INTRODUCTION 

Hong-Ou-Mandel interference is a two-photon quan- 
tum interference effect whereby two indistinguishable 
photons incident at either port of a balanced beam spht- 
ter will always be detected to coalesce at one or the 
other output port of the beam splitter, but never with 
one photon at each output port. This interference ef- 
fect diminishes as the distinguishability of the photons is 
increased. It was first demonstrated by Hong, Ou and 
Mandel in 1987 using photons produced via parametric 
down-conversion |2] . The importance of this effect in op- 
tical implementations of quantum information process- 
ing schemes was realized later [31 U , motivating further 
demonstrations. They have been achieved both by us- 
ing a single source [H [6] and by using two independent 
sources [7HS] of photons. In addition to these pulsed in- 
terference experiments, two-photon quantum interference 
using continuously-driven sources of resonant fluorescent 
light has also been demonstrated [10]. 

However, until the recent experiment of Lang et al. [1] , 
a comparable demonstration of the same effect with 
microwave-frequency photons had been lacking. This 
experiment was performed using a circuit QED system 
|lll I12j , configured as pulsed sources of single microwave 
photons [T3j. Figs. 2(c)-(f) of Ref . [T| show both the mea- 
sured and calculated second-order correlation functions 
of the microwave fields at the output of the beam splitter. 
The calculation of these theory curves, corresponding to 
inputs consisting of trains of Lorentzian photons, is pre- 
sented here, along with the same calculation for Gaus- 
sian photons. One could also perform the experiment 
with continuously-driven circuit QED systems, operated 
in the regime of resonant photon blockade [14] . We also 



calculate the output second-order correlation functions 
in this case, though these are not directly relevant to the 
experiments described in Ref. [TJ 



The calculation of the correlation functions in the 
pulsed case proceeds in a manner similar to that of Legero 
et al. |15j . undertaken to describe their experiments on 
time-resolved two-photon quantum interference [6]. We 
generalise their calculations in a number of ways. First, 
we consider Lorentzian (in addition to Gaussian) single 
photons. This is because here the spatio-temporal prop- 
erties of the emitted photons are determined purely by 
the decay rate of the circuit QED source, rather than by 
additional broadening mechanisms. Second, we calculate 
the correlation functions for a train of single photons, 
leading to features in the calculated correlation functions 
at (and around) integer multiples of the pulse period. 
Further, these correlation functions are evaluated as cir- 
cular correlation functions, meaning that the correlations 
for a finite pulse train are obtained by wrapping the time 
index around the finite sequence of pulses. Third, in or- 
der to make a comparison with experimental data, we 
subject the calculated correlation functions to filtering 
representative of the filtering that the measured data is 
subjected to. 



The second-order correlation functions at the output 
of the beam splitter are evaluated, in terms of correlation 
functions of input fields, in Sec. IH] These are explicitly 
evaluated, for both individual and trains of pulsed single- 
photons (Gaussian and Lorentzian), in Sec. HI The cor- 
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relation functions for continuously-driven sources, both 
without and with polarisation degrees of freedom, are 
evaluated in Sec. II VI Note that neither the details of the 
extraction of the microwave correlation functions from 
the directly measured correlation functions in Ref. T, nor 
the details of the quantum state tomography employed 
in Ref. [U are described here. They have been discussed 
elsewhere [TOIITT] . 



II. OUTPUT SECOND-ORDER CORRELATION 
FUNCTIONS 

A. Intensity Correlation EYinctions 

Suppose that the two fields incident on the input ports 
of the beam splitter are described by dimensionless, time- 
dependent quantum fields denoted by a'{t) and b'{t). 
Unitarity of the (balanced) beam splitter implies that 
the output fields, a(t) and b{t), are then given by the 
Heisenberg picture evolution [18], 



ait) 

m 



1 

72 



1 -1 
1 1 



a'(t) 
b'it) 



(1) 



We are interested, primarily, in the statistics of the pho- 
tons output from the beam splitter. These are described 
by correlation functions of the number operators of the 
output modes, the number operators being denoted by 



flcit) 



'(t)c(t). Then, the quantities we seek are 



the normally-ordered correlation functions of the output 
mode number operators; 

(: fiait + T)hr{t) :) = {f\t)a\t + T)a{t + T)f{t)) , (2) 



with f representing either a (for the output intensity 
auto- correlation function) or b (for the output intensity 
cross-correlation function). These correlation functions 
allow one to discriminate between quantum and classical 
fields [18] . Since the intensity of a field is proportional to 
the number of photons it contains, henceforth we shall 
refer to correlation functions of the form of Eq. (pi) as 
intensity correlation functions. The term second-order 
correlation Junction shall be reserved here for functions 
of the difference in detection times (t) alone, to be intro- 
duced below. The correlation functions (: n(,(t-|-r)nr(i) :) 
give the same results due to the symmetry of Ec^. (flj) . 

Clearly, we can obtain the intensity correlation func- 
tions of Eq. ^ in terms of the input fields by substi- 
tuting the expressions in Eq. (fTl). The result includes 
sixteen fourth-order moments of the time-dependent in- 
put fields. Now, if it is assumed that the sources pro- 
duce separable (i.e. not entangled) fields (clearly a reason- 
able assumption for independent sources), the four-term 
cross-correlation functions of input fields may each be 
decomposed into products of two auto-correlation func- 
tions of the input fields. Therefore, the output intensity 
auto-correlation and cross-correlation functions are: 
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(3) 



where the + and — signs correspond to the auto- 
correlation and cross-correlation functions, respectively. 
Note that the summation is over the set of permuta- 
tions of the input mode operators a' and b'; that is, 
n(a',6') = {{a' ,b'),{b' ,0,')}. The output intensity corre- 
lation functions depend on the input intensity correlation 
functions [first term in Eq. ([3|)], the products of the inten- 
sities of the two inputs (second term) , the products of the 
first-order correlation functions of the two inputs (third 
term), and on a number of phase-dependent correlation 
fmictions of the inputs (terms four to eight). 



B. Phase-Dependent Moments 

For the purpose of demonstrating Hong-Ou-Mandel in- 
terference, the contributions of the fourth through to 
the eighth terms to the intensity correlation functions of 
Eq. ([3|) are neither necessary nor desirable. We will con- 



sider their contribution using the example of the fourth 
term, though the same arguments may be applied to the 
fifth through to the eighth terms. Suppose that the fields 
corresponding to the input modes, c and d in Ec^. (|3|, 
can be written as slowly-varying fields modulating os- 
cillations at carrier frequencies centred on wg and with 
a difference A, and with some relative phase offset 6: 
c{t) = c(t)e+*("°+'^/2)t+»e ^^d J(^) ^ j(^)g+»(^„-A/2)t^ 

Accordingly, we can rewrite the fourth term in Eq. ([3|) as 



{C^ [t + T)cHt)) {d{t + T)d(t)) 
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If A is non-zero this expression will remain oscillatory 
even in the limit t — >■ oo. However, if |A| is large com- 
pared with the inverse time-scales of the features of in- 
terest in the input fields, this expression will average to 
zero over the time-scales of interest. Alternatively, for 
single-photon input states the phase-dependent second 
moments will always evaluate to zero. Even if the expres- 



sion in Eq. Q does not automatically vanish, its effect 
can be averaged to zero over multiple runs of the exper- 
iment if the relative phase 9 is uniformly distributed. In 
all cases of interest to us, the phase-dependent moments 
will vanish and the expression in Eq. (l3| reduces to: 

(: na{t + T)ha(b){i) ■) 
= \ E [{^{t)cHt + T)c{t + r)c{t)) 

(c,d)en(a',b') 

+ {c\t)c{t)){dHt + T)d{t + T)) 

± {c\t + T)c{t)) {d\t)d{t + r))]. (5) 

Here the bar over the correlation function indicates that 
it is to be evaluated over many runs of the experiment. 
Henceforth, this shall be assumed for the correlation 
functions quoted and we drop this notation. The output 
intensity correlation functions are completely determined 
by the input intensity correlation functions and products 
of the input first-order correlation functions. For the out- 
put intensity cross-correlation function, the two-photon 
interference is seen in the cancellation of the second and 
third terms at r = in Eq. ([s]). On the other hand, the 
second and third terms will add constructively for the 
output auto- correlation function at r = 0, corresponding 
to photon coalescence. Whether or not this constitutes 
Hong-Ou-Mandel interference in the usual sense is deter- 
mined by the input intensity correlation functions in the 
first term of Eq. (Is]) . Although the preceding statements 
are true in general, the signature of the interference for 
non-zero r will depend on the form of the input fields. 



C. Second-Order Correlation Functions 

The output intensity correlation functions, as ex- 
pressed in Eqs. ([3]) and ([5]), are explicitly dependent on 
two detection times, t and t-\-T. The second-order corre- 
lation functions, giving the likelihood of detecting a pho- 
ton at some time r assuming detection at time zero, may 
be obtained either by integrating or by taking a limit: 



.(2)/ 



Gcdi^) = J dt{:n,{t + T)nd{t):) or 



(6) 



respectively. Whether one integrates or takes a limit de- 
pends on whether the system's evolution is described in a 
(non-rotating) laboratory or a rotating frame. The quan- 
tities in Eq. ([6]) are conventionally referred to as second- 
order correlation functions |18j . Similarly, we introduce 
the first-order correlation functions: 



(7) 



G«M = j dt{cHt)c{t + T)) 
G^P{t) = \\m{c\t)c{t + T)). 



In the case where the second-order correlation functions 
are obtained by taking limits, the output second-order 
correlation functions follow immediately from Eq. (p]) as: 
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where we have introduced the short-hand notation 



.(2)/ 



-(2) 



Gc (t) = Gcc (t). If both sources exhibit identical 
statistics, meaning GI^/{t) = Gi,'{t) = G' '{t) and 

G^^\t) = G^^\t) = G['^\t), the output second-order 
correlation functions are simply: 
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As for Eq. ([5]), the second and third terms of the second- 
order cross-correlation functions cancel at r = 0, corre- 
sponding to photon interference at the beam splitter. We 
now consider the evaluation of these correlation functions 
for particular forms of the input fields. 



III. PULSED SINGLE-PHOTON SOURCES 

A. Spatio-temporal Modes and Correlation 
Functions 

The case relevant to the experiment described in Ref. [T] 
is that in which each source produces a highly-pure, 
microwave-frequency single-photon state. We shall eval- 
uate the correlation functions of Eq. (p]) by introducing 
spatio-temporal mode functions to describe these single 
photons. First, the time-dependent quantum fields intro- 
duced above, c(i) in general, may be represented in the 
frequency domain as c(w) = -^ J dt e+"^*c(t). Then 
a single-photon pulse is conventionally represented using 
the spectral density <&(w); that is, 



1) = duj ^{Lo)c\io)\0) 



(10) 



One can show, by considering the action of c(t) on the 
state |1) , that the time-dependent quantum field c{t) may 
be represented as 



c{t) = ({t)co, where C(^) 



1 



/2n 



duj ^{Lo)e 



(11) 

is the spatio-temporal mode function of the single-photon 
pulse, and cq is the annihilation operator correspond- 
ing to this spatio-temporal mode. It is then clear that 
(c'{t)c{t)) = |^(t)| , giving the probability of the photon 
being detected at time t. 

In the experiment of Ref. 1, there are two indepen- 
dent microwave-frequency sources producing highly-pure 



single-photon states. Accordingly, we can write the beam 
splitter input fields (from each source) in terms of spatio- 
temporal mode functions as 



a\t) = Ca{tK, y{t)^Cbm- 



(12) 



Accordingly, the two-mode state input to the beam split- 
ter may be written as d''''(t)6''''(t) |0). From Eq. ([5|, the 
output intensity correlation functions are: 



(: 7ia(< + T)na(f,)(^) :) 

= ]\Ca{t + T)Cb{t)±Ct{t + T)Ca{t)f 



(13) 



Here, as in Eq. (5|, the potential for cancellation of the 
terms on the right-hand-side in the case of the intensity 
cross-correlation function is clear. Indeed, at r = 0, we 
see that the intensity cross-correlation function vanishes 
irrespective of the precise form of the spatio-temporal 
mode functions. Since there are now no additional con- 
tributions to the output intensity cross-correlation func- 
tions, the intensity correlation function vanishes at r = 0. 
This suppression is the signature of Hong-Ou-Mandcl in- 
terference. On the other hand, the terms in Eq. ( 13 ) add 



constructively for the intensity auto-correlation function, 
indicative of photon coalescence. Of course, the form of 
the intensity correlation functions for non-zero r depends 
on the form of the mode functions, and we evaluate this 
now for Gaussian and Lorentzian photons. 



B. Single Gaussian Photons 

One form of mode function commonly produced by 
single-photon sources are Gaussian modes. For exam- 
ple, this is the case for trapped atom sources in which 
linewidth broadening mechanisms are dominant |B]. In 
the subsequent section we consider the case of Lorentzian 
mode functions, which arise naturally for circuit QED 
systems. Given the recently demonstrated ability to 
shape the waveforms of single photons [H] , we could con- 
sider more general mode functions as well. Two Gaussian 
modes, with a carrier frequency difference A (centred on 
Wo), a relative temporal offset St (centred on t — 0), and 
a common pulse width specified by a, are described by 
the spatio-temporal mode functions: 



C.wW = V^ 



cxp 



{t±ST/2Y 

2^^2 



-i{ujo±A/2)t 



(14) 
Substituting Eq. ((T4| into Eq. ( 13 ) gives the explicit 



form of the output intensity correlation functions. Sub- 
sequently, integrating over the first detection time t leads 
to the second-order correlation functions, 

(2) ^ cosh (r (5t/ct^) ± cos (rA) 

aaiab)^^) - 2V2TTa^ 



The total correlation probabilities^ Paa and Pab, are 
then obtained by integrating over the detection time r. 
The correlation probabilities Pcd simply give the proba- 
bility of detecting one photon in mode c and one photon 
in mode d, for one pair of input photons, irrespective of 
the detection times. They are: 
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aai^ab) 



1 ± exp 



2^2 



A^a 



5t^ 
2^2 



(16) 



Note that Pba = Pah and Pbb = Paa, such that 
J2c de{a b} ^cd = 1, as required for a beam splitter that is 
assumed to be lossless. The total probability for photon 
coalescence at one output port is Paa + Pbb = '^Paa , while 
the total probability for one photon in either output port, 
the coincidence probability, is P^ = Pab+Pba = "^Pab- The 
correlation probabilities may be regarded as probability 
distributions over the controlled parameters 5t and A. 
For indistinguishable photons {5t = and A = 0) the 
coincidence probability is Pc = 0, and for fully distin- 
guishable photons (i.e. \5t\ ^ a and/or |A| :^ 1/cr) we 
have Pc — >■ 1/2. The coincidence probability, plotted as 
a function of either 5t or A, exhibits the classic Hong- 
Ou-Mandel dip j2j. 



C. Single Lorentzian Photons 

The same calculation may be performed for Lorentzian 
photons. In the experiment of Ref. |Ijthe spatio-temporal 
mode function of the released photon is determined by 
the decay rate of the coupled circuit QED system. Ac- 
cordingly, these photons have the shape of a truncated 
exponential decay. Two such mode functions, with a rel- 
ative temporal offset 5t and a carrier frequency difference 
A, are given by: 



Caib)it) = y2^e-^(*±^*/2)w(t±5i/2)e- 



i{u)o±/\/2)t 



(17) 



where 7 is the field decay rate of the coupled system, wg 
is the frequency about which the two carrier frequencies 
are centred, and u{. . .) denotes the unit step function. 
Note that the field decay rate is the characteristic decay 
rate of the electric field, in contrast to the commonly 
quoted (cavity) decay rate of the energy stored in the 
cavity. According to Eq. (171, the first photon (in mode 
a') is "released" at t = —dt/2, while the second photon 
(in mode b') is released at t = +6t/2. In the frequency 
domain, the mode functions correspond to the spectral 
densities, 



'^a{b){uj) 
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TT 7 - i [w - (cjQ ± A/2)] 
X exp {-i [uj - {uo ± A/2)] dt/2} . (18) 



This representation makes it clear why the exponentially- 



exp [— ((5t^ + T'^) /2(T^] . (15) decaying single-photon pulses described by Eq. (17) are 



referred to as Lorentzian photons. Of course, the mode 
function thus defined is unphysically sharp. The filter- 
ing effects of a finite measurement bandwidth shall be 
addressed below. 

Substituting Eq. (Tf^ into Eq. ( 13 ), assuming that 6t > 
and r > 0, we find the output intensity correlation 
functions to be: 



(: fiait- 
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(19) 



where the upper (lower) result corresponds to the auto- 
correlation (cross-correlation) function. Eq. ( 19 1 consists 
of two contributions, corresponding to different combi- 
nations of the unit step functions switching on. The 
first product of step functions in each switches on when 
both the first detection (i) is after the release of the first 
photon {—St/2) and the second detection time (i -I- t) 
is after the release of the second photon {+St/2). The 
second unit step function is switched on when the first, 
and therefore the second detection time (t > 0) is after 
the release of the second photon. Integrating over the 
first detection time t leads to the output second-order 
correlation functions. 
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(20) 



Again, Eq. (20 1 is split into two terms, depending 



on whether the difference in measurement times (r) is 
greater or less than the temporal delay (St) between the 
release of photons from either source. Integrating over 
the difference in detection times r (positive and nega- 



tive) we obtain the correlation probabilities, 
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As before, we have Pi 



bb 



J^^n.n 1 J^^hi 



Pab, EPcd - 1, 

and the coincidence probability is Pc = Paa + Pbb — '^Paa- 
Again, these probabilities may be regarded as probability 
distributions over St and A. If the photons are indistin- 
guishable {St = and A = 0), Pc = 0, while for fully 
distinguishable photons {St 3> I/7 and/or |A| ^ 7), 
Pc ^ 1/2. 

D. Trains of Single Photons 

The experiment in Ref. Jl is performed with trains of 
pulsed single photons. In general, this can lead to fea- 
tures in the correlation functions at (and around) inte- 
ger multiples of the pulse period. These are evaluated 
here. In order to determine the output intensity cor- 
relation functions in this case, we consider each photon 
as belonging to an independent spatio-temporal mode, 
and assume that the state of the many-photon system 
is separable with respect to this basis of spatio-temporal 
modes. Accordingly, the time-dependent quantum fields 
at the inputs may be decomposed as the sum of the fields 
corresponding to each spatio-temporal mode: 
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with the spatio-temporal mode operator commutation re- 
lations being [cq^., Cqj,,] = Si-^', where c' denotes an input 
mode lowering operator, a' or b' . The appropriate spatio- 
temporal mode functions for the k photon of each input 
(assumed Lorentzian here) are: 

Xg-7(t-fet^±5t/2)y (t - ktp ± St/2) ,(23) 

with tp denoting the pulse period. Substituting the field 
decompositions of Eq. ( 22 ) into the intensity correlation 
functions of Eq. ^ leads to: 
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FIG. 1: Output second-order (a) cross-correlation functions and (b) auto-correlation functions for trains of Lorentzian single 
photons input at each port of a beam splitter with pulse period tp. The delay of one source with respect to the other in each 
case is 5i = ns (black line), 50 ns (dashed blue line), 100 ns (dashed red line) for tp — 0.5 /is. For indistinguishable photons 
(black line), the signal in the first pulse period of the cross-correlation function is suppressed, indicating interference. 



Using the assumption that the input state is separable 
with respect to the basis of spatio-temporal m.odes 



introduced, Eq. ( 24 1 reduces to 
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(25) 



The output 
G 



(2) 



second-order correlation 
then obtained 



functions, malised such that the peaks at ±ktp in GJ^j, (r) are unity 



aafah)(^)' ^^^ then obtained by integrating over 
the first detection time t, as in Eq. ^. 

These correlation functions are plotted in Fig. [T] here 
and Figs. 2(c)-(f) of Ref. fTl Note that we actually cal- 
culate circular correlation functions, meaning that we 
evaluate these sums for a finite pulse train, and then 
wrap the time index around the pulse train. This means 
that the peak heights at later pulse periods will be at the 
same height as those at earlier pulse periods, rather than 
showing an artificial decay in peak height that would oth- 
erwise be produced due to the finite nature of the pulse 
train. Furthermore, the correlation functions are nor- 



for indistinguishable photons input. 

As expected, the signature of interference in Fig. [l] 
is the suppression of the output second-order cross- 
correlation function for small delays r for indistinguish- 
able photons {St = 0). As the distinguishability of the 
photons is increased through increasing St, one clearly 
sees corresponding peaks in the first repeat period. The 
peaks in the correlation functions at ±ktp {k an integer) 
exhibit side peaks at ±5t. This is easily understood: a 
given photon of each pair may be correlated will be cor- 
related with the first or second photon of the subsequent 
pair, as well as with the first or second photon of the 
preceding pair. 



E. Filtered Response 

Integrating the intensity correlation function of 
Eq. (25) results in second-order correlation functions 
with the very sharp features seen in Fig. [T] due to the 
abrupt switching of the single-photon pulses described 
by the unit step and exponential function in Eq. (17 1. 



However, such features are not observed in experimental 
data due to filtering in the detection scheme. To facilitate 
comparison with experimental data, we incorporate this 
filtering into our calculations. The filtering is performed 
on the time-dependent signals obtained by the data ac- 
quisition scheme in Ref. IT] denoted by Sa(b){t)- The 
second-order correlation functions of the fields are ob- 
tained from the filtered version of these signals, S^fj^Jt), 
given by the convolution integral. 






w = 



+ CX3 



Saib){t')f{t-t')dt', 



(26) 



where / is a filter function in the time domain. Now the 
correlation functions of the output fields of interest may 
be obtained from the correlation functions of these mea- 
sured signals, and so for the sake of calculations, we may 
implement the filtering directly on the spatio-temporal 
mode functions; that is, 



+ 00 



Cmkii)= Caib){t')f{t-t')dt'. 



(27) 



The filtering in Ref. l!was implemented in two stages. In 
the first stage, a square window filter was employed to 
eradicate specific frequency components associated with 
the analog-to-digital conversion. A subsequent filtering 
stage, implemented as a finite-impulse-response Cheby- 
shev filter [20|, reduced other noise sources. For the pur- 



pose of calculations, Eq. (27) was implemented only once. 



but with a filter function representative of that which was 
done experimentally. The intensity correlation functions 
are then given by Eq. (25), evaluated using the filtered 



spatio-temporal mode functions of the form of Eq. (27) 



The second-order correlation functions are then obtained 
by integrating over the first detection time t, as in Eq. (pi), 
and are shown as the theory curves in Figs. 2(c)-(f) of 
Ref.ffl 



IV. CONTINUOUSLY-DRIVEN SOURCES 

A. Introduction 

Rather than assuming the circuit QED systems are 
configured as pulsed single-photon sources, we can as- 
sume that they are continuously-driven and operated in 
the regime of resonant photon blockade [14]. The out- 
put microwave field from each source will exhibit sub- 
Poissonian and anti-bunched photon statistics [HJ E2] . 
leading to the possibility of observing Hong-Ou-Mandel 



interference at the beam splitter [TU]. Treating only the 
lowest transition of the coupled circuit QED system [11] , 
we consider the source simply as a driven quantum two- 
level system (TLS). The physics of each source is then 
essentially the physics of resonance fluorescence [531 IMl j 
with additional dephasing included. Note that this TLS 
is in fact formed from the coupled system, and is therefore 
not the "qubit" typically employed for superconducting 
quantum information processing. 

We denote the TLS transition frequency by Wc (where 
c denotes the corresponding input mode, either a' or 
&'), and assume that it is subject to driving by a res- 
onant monochromatic, large-amplitude microwave field. 
Each system is described, in a frame rotating with re- 
spect to Hq = [hiOc/2)az, by the Hamiltonian Hj = 
h (51a+ -f 51*cr_) /2 where fl = 2g/3 is the Rabi frequency, 
with /3 denoting the coherent amplitude of the driving 
field and g being the Jaynes-Cummings coupling between 
the TLS and the driving mode [llj. The coupling of the 
TLS to its environment may be fully described, in the 
Markovian white noise approximation, by damping at a 
rate 71 and dephasing at a rate 7^ . Accordingly, the evo- 
lution of the TLS density operator p may be given in the 
form of a Lindblad master equation [18] as: 



[na+ + n*(7^,p] + 712? [cr_]/? 



7p 



(28) 
where V\A\ = ApA^ - \A^ Ap - \pA^ A is the so-called 
dissipative superoperator. Note that in this model each 
source is fully characterized by the parameters 7i,7p, SI, 
and LUr- 



B. Output Correlation Functions 

As in the case of pulsed sources, the quantities of in- 
terest are the second-order correlation functions of the 
microwave fields at the outputs of the beam splitter. The 
master equation of Eq. ( 28 ) facilitates the evaluation of 



correlation functions of TLS operators, and from the the- 
ory of atomic spontaneous emission [24] . we know that 
the output field will be proportional to the TLS lowering 
operator. 



m 



Xt), 



(29) 



where c denotes the input mode a' or b' . The output 
intensity correlation functions then follow from Eq. (|3|, 
to within a constant factor, by making the substitution 
in Eq. ( 29 ) . Rather than explicitly evaluating the coeffi- 



cient in Eq. ( 29 ), we shall ultimately calculate normalised 



(field) correlation functions for which the coefficients will 
cancel. 

First we calculate the output correlation functions in 
terms of the input correlation functions. Suppose that 
the TLS transition frequencies are uiq ± A/2 (for the 
sources corresponding to the input modes a' and 6', re- 
spectively), and that we are working in an interaction 



picture with respect to Hq — ^^f^ta' b'} (^o/2)cr^- Then, 
from Eqs. (Is]) and (29), the output intensity correlation 
functions (to within a constant factor) are: 

(: 7ia(t + r)7i„(t)(i) :) 

-I E [K(iK(t + r)a^(i + T)ai(t)) 

(cM)en{a' M') 

+ «(i)a^(i))(4(i + r)al(t + r))] 

±e-^(a^ (t + T^: (t)) «' (t)a^' (t + r)) 

-"^K'(t + r)a^'(t))(4'(t)a^'(t + r))". 

(30) 



1 r 
4 



Note that the vahdity of this expression depends on 
averaging out exphcitly phase-dependent contributions 
by randomising the relative drive phases over repeated 
runs of the experiment, as discussed in Sec. |IIB| Taking 
the limit i — >• oo of Eq. ( 30 1 we have the output second- 



order correlation functions: 



G„ 



iab)ir) - ^{Gi^^(r) + Gi?)(r) + 2GiV(0)GW(0) 



±e-^Gi!'(r)[GW(r)' 



±e 



-itA 



GiVM]*G«M}, (31) 



where the correlation functions on the right-hand-side are 
now TLS correlation functions, 



GfKr) 



yiuY{a%{t)a'L{t + T)), 



lim {(T%[t)<j%{t + T)<j1{t + T)a'L{t)). 

t—^oo 



(32b) 



If one assumes that (carrier frequency aside) the sources 

(2'] (2) 

have identical characteristics, meaning G\, — G^, = 



Cf^ and G^} 



-(2) 



G 



(1) 



Gj- , then we have: 



G", At) 

aa[ao) ^ ^ 



iJGf)(r)±cosrA|GW(r)" 



Gf)(0) 



(33) 



Again, the second and third terms in the output cross- 
correlation function of Eq. (33) cancel at t = 0. For non- 



zero r, the interference is oscillatory, an effect sometimes 
referred to as a quantum heat [3]. 



C. Resonant Photon Blockaded Inputs 

Now we turn to the evaluation of the input TLS cor- 
relation functions. The system of equations for the ex- 
pectations of TLS operators are readily calculated from 



Eq. ( 28 ) . They are commonly referred to as the optical 



Block equations^ and may be expressed as the inhomoge- 
neous system [18], 



dt 



(a) ^A{a)+ b, 



(34) 



where we have introduced the vector of expectations of 
TLS (Pauli) operators, {a) = {{a+) , (a^) ,{a^)f . The 
inhomogeneity is given by 5 = (0, 0, —71) , and the sys- 
tem matrix is given by 



-72/2 -fiA -in*/2 

-72/2 - iA in/2 

—iQ ifl* —71 



(35) 



where we have also introduced the total dephasing rate, 
72 = 7i + 27V One can easily find the steady-state solu- 
tion of Eq. (34), (CT)gg, and subsequently recast Eq. (34) 



as the homogeneous system, d (a) /dt = A (a) where 
((?)' = ((?) — ((?)^^. Now the time-dependent solution of 
the homogeneous system may be obtained by diagonal- 
ization, and subsequently the solution to Eq. (34) may 
be written in the form, 



(a(i + r)>-G(r)(a(i))+d(r), 



where we have introduced the notation 



X-^AXt v-l 



d[{T) = (l-Xe^"^^-)(a),„ 



(36) 

(37a) 
(37b) 



(32a) with X being the matrix of eigenvectors of A 



The form of the solution in Eq. ( 36 ) is amenable to 



application of the quantum regression theorem |24| . In 
particular, for an clement of the solution in Eq. (36), we 
can write 

(a,(t + r)> =EG.,(T)(a,(t))-Kd.(T), (38) 



where i G {+, — , z}, the summation index j runs over this 
set, and Gij{T) and di{T) denote elements of the matrix 
G{t) and vector d{T), respectively. The quantum regres- 
sion theorem allows us to evaluate correlation functions 
as: 

{ak{t)a^{t + T)aiit)) = ^ ^u W('^fcW^j W^/W) 

j 
+d,{T){<Jk{t)<Ji{t)}. (39) 

The TLS correlation functions in Eq. (|30| are obtained 



using Eq. ( 39 ) and the algebraic properties of Pauli op- 
erators. 



According to Eq. ( |33[ ), the output second-order corre- 
lation functions are fully determined by the input first- 
order and second-order correlation functions. We can 
write out G^ (r) and G^ (t), which are well-known re- 
sults from the physics of resonance fluorescence |Mj . For 




^ aa V ' J 




0.5 







If 



A/27r = 

A/27r = 31.4 MHz 

A/2n = 62.8 MHz 
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FIG. 2: Output second-order (a) cross-correlation functions and (b) auto-correlation functions for each source continuously- 
driven in the regime of resonant photon blockade. Again, interference is observed in the suppression of the cross-correlation 
function for small r. Photon distinguishability may be introduced through a carrier frequency difference A, which leads to an 
oscillation in the observed interference. 



the first-order correlation functions we have: 

|2 



G['\r) = 



7172 + 2 \n\ 
1 



7? 



-72T/2 



7i72 + 2|f^| 



(27i+72)r/4 



e -"' '-■'-(A+e+'^^ + A-e-''^) 



(40) 



where we have introduced the notation 

4k = y(2^^^^^Jf^T(^, (41a) 

A± = {2|r!p-272+7i72±[2|f]|2(67i-72) 



71 (271 - 72)^ /4k| / (7172 + 



2|r2| 



(41b) 



Note that the usual resonance fluorescence spectrum fol- 
lows from the result in Eq. ( 40 ) , and in the strong driving 



regime {\fl\ ^ 71,72) we can obtain a simple analytical 
expression describing the so-called Mollow triplet [531 HI] 
feature in the spectrum. For the second-order correlation 
function we have: 



Gf'ir) = 



m 



7i72 + 2|f^| 



I _ g-(27l+72)T/4 



X ( cosh KT H sinh kt 

4k 



.(42) 



It is more common to quote the normalised result, which 

is: 



Gf'(O) 



= 1 _ p-(27l+72)T/4 



= 1 



2ti ~I~ ^9 

X ( cosh KT -\ sinh kt 1 . (43) 

4k ' 



Of course, this is a normalised TLS correlation function. 
Here we shall normalise the output field correlation func- 
tions by imposing the requirement that they asymptote 
to unity. 

Now the output second-order correlation functions are 
given by Eq. (33), using the results in Eqs. (40) and (42 1. 
They are shown in Fig. [2J As expected, the interfer- 
ence effect is observed in the suppression of the cross- 
correlation function for small delays r, while introduc- 
ing distinguishability through the frequency difference A 
leads to an oscillation in the observed interference for 



non-zero r. 



D. Photon Polarisation 

In the continuously-driven case just considered, pho- 
ton distinguishability can only be introduced through a 
difference in the carrier frequencies. An alternative ap- 
proach, relevant in a three-dimensional circuit QED ar- 
chitecture |25| . is to introduce distinguishability through 
a polarisation degree of freedom. To treat this problem, 
we decompose the time-dependent fields into their (hor- 
izontal and vertical) polarisation components: 



Z{t) ^ Ch{t)+C^{t), 



(44) 



where the subscripts h and v denote the horizontally and 
vertically polarised components, respectively, of the field. 
Generalizing the expression of Eq. (pi) to include polari- 
sation indices, and then summing over these indices, we 
can write the output intensity correlation functions as: 

{: na{t + T)fi^{t) :) 
= E (4Wa|(i + T)a,(t + r)cfe(t)), (45) 

k,l^{h.,v} 

with c representing either output mode (a or 6), as before. 
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FIG. 3: Output second-order (a) cross-correlation and (b) auto-correlation functions for each source continuously-driven in 
the regime of resonant photon blockade, for a range of polarisations of one input with respect to the other (<^ = 0°, 45°, 90°). 
Here the carrier frequencies of both inputs are assumed to be the same. For indistinguishable photons (0 = 0°), the interference 
signature of a suppressed cross-correlation function for small r is observed. As the indistinguishability (i.e. relative polarisation) 
is increased, the interference effect is reduced. 



Now we take the components of the input fields to be 
related to the TLS lowering operators by: 






al{t), aUi) = 0, 



(46a) 



a''^ (t) cos (j), b'^ (i) - cr^' (t) sin 0. (46b) 



That is, the polarisation of one input with respect to 
the other is specified by the angle (j). Accordingly, using 
Eq. ^ and the decomposition of Eq. ( 44 1 , we can find 
the polarisation components of the output fields. Sub- 
stituting these into the generalised intensity correlation 
functions of Eq. (|45| we find: 



(: fiait + T)na[b){t) :) 

(c,d)en(a',6') 

+ {a%{t)at{t)){o<i{t + T)a'L{t + T)) 

± cos" 4>{(rl{t + T)at{t)) {ai{t)aUt + r))} .(47) 

The photon distinguishability, and therefore the presence 
or absence of two-photon quantum interference, is deter- 
mined by the relative polarisation angle (f>. For = 0, 
the photons at either input are indistinguishable, and the 
second and third terms in Eq. (47) cancel for the cross- 



correlation function, corresponding to interference. For 
the auto-correlation function, the second and third terms 
will add constructively, indicative of photon coalescence. 
For (p = 7r/2, the photons at either input are fully distin- 
guishable, and the last term in Eq. (47) makes no con- 
tribution. Now taking the limit i — > cxd, as per Eq. (p|, 
yields the second-order correlation functions, 



-(2) 



aa[ab) V / 



1 

4 



E {G(^)(r) + GW(0)G«(0) 

n(a',fc') 

.2 0GW(r)[GW(r)]*}. (48) 



(c,d)en(a',fc') 

±cos 



Under the additional assumption that the two sources 
(polarisation aside) have identical properties, meaning 

G^'\r) = G<;^\r) ^ G^^\r) and Gi]\r) = G^^\t) ^ 
.(1) 



G\ (t), then we find: 

Gaa{ah){T) ^ 



^(2) 



Gr\T)+ G^i'iO) 



.(1)/ 



± COS^ ( 



G^'ir) 



(49) 



As for Eq. (47), the dependence of photon interference 
on the relative polarisation angle 4> is clear. These out- 
put correlation functions are shown in Fig. |3] Clearly, 
access to a polarisation degree of freedom allows one to 
smoothly interpolate between the cases of distinguishable 
and indistinguishable photons in the continuously-driven 
case. 



V. CONCLUSIONS 

Two-photon quantum interference at a beam splitter 
is observed in the suppression of the output second-order 
cross-correlation function for small detection delay times. 
In the case of pulsed single-photon sources, this quan- 
tity may be integrated over all detection delay times to 
reproduce the classic Hong-Ou-Mandel dip in the co- 
incidence probability at the output. Accounting for a 
train of photons input leads to sidebands on the second- 
order correlation functions around integer multiples of 
the pulse period, while the inclusion of a finite detection 
bandwidth broadens the sharp features in the correlation 
functions that would otherwise be observed. The second- 
order correlation functions calculated are presented in 
Figs. 2(c)-(f) of Ref. [U showing very good agreement 
with the experimental data. In the case of continuously- 
driven sources, photon distinguishability may be intro- 
duced through a difference in carrier frequency or via a 



polarisation degree of freedom. In the former case, an 
oscillation in the quantum interference is expected. 
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